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ABSTRACT 



Disks around young stars are known to evolve from optically thick, gas-dominated protoplanetary disks to optically thin, almost gas- 
free debris disks. It is thought that the primordial gas is largely removed at ages of ~ 10 Myr and indeed, only little amounts of gas 
have been deduced from observations for debris disks at ages of S 10 Myr However, gas detections are difficult and often indirect, 
not allowing one to discern the true gas densities. This suggests using dynamical arguments: it has been argued that gas, if present 
with higher densities, would lead to flatter radial profiles of the dust density and brightness than those actually observed. In this paper, 
we systematically study the influence of gas on the radial profiles of brightness. We assume that dust is replenished by planetesimals 
orbiting in a "birth ring" and model the dust distribution and scattered-light brightness profile in the outer part of the disk exterior to 
the birth ring, under different assumptions about the gas component. Our numerical simulations, supported with an analytic model, 
show that the radial profile of dust density and the surface brightness are surprisingly insensitive to variation of the parameters of a 
central star, location of the dust-producing planetesimal belt, dustiness of the disk and — most importantly — the parameters of the 
ambient gas. The radial brightness slopes in the outer disks are all typically in the range -3. ..-4. This result holds for a wide range of 
gas densities (three orders of magnitude), for different radial profiles of the gas temperature, both for gas of solar composition and gas 
of strongly non-solar composition. The slopes of -3. ..-4 we find are the same that were theoretically found for gas-free debris disks, 
and they are the same as actually retrieved from observations of many debris disks. Our specific results for three young (10-30 Myr 
old), spatially resolved, edge-on debris disks (yS Pic, HD 32297, and AU Mic) show that the observed radial profiles of the surface 
brightness do not pose any stringent constraints on the gas component of the disk. We cannot exclude that outer parts of the systems 
may have retained substantial amounts of primordial gas which is not evident in the gas observations (e.g. as much as 50 Earth masses 
for p Pic). However, the possibility that gas, most likely secondary, is only present in little to moderate amounts, as deduced from gas 
detections (e.g. ~ 0.05 Earth masses in the /? Pic disk or even less), remains open, too. 
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1. Introduction 

In the course of their evolution, circumstellar disks transform 
from optically thick, gas-dominated protoplanetary disks to op- 
tically thin, almost gas-free debris disks. How the gas is re- 
moved is not known in detail, but the removal is thought to be 
the result of a quick inside-out process at an age of ~ 10 Myr 
jXTexander 2008; Hillenbrand 2008). It may be either due to the 
UV switch mechanism resulting from an interplay between pho- 



200C; 



toevaporation and viscous accretion (e.g.'Hollenbach et al 
Clarke et al. 2001; Takeuchi & Lin 2005; Takeuchi et af 
Alexander & Armitagel I2007D or due to gap opening by hid 
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den g iant planets (e.g. iLubow et"aLlll999l; iLubow & D'Angelol 
120061) . Although the former effe ct seems to be slig htly preferred 
by observational statistics (e.g. ICieza et aL I l2008h . it is not yet 
possible to distinguish between thern with certainty (iNaiita et al.l 
I2007t lAlexaiidea2008t lHrilenbrandl2008h . 

In this paper, we concentrate on a more advanced phase of 
system's evolution: the debris disk stage. Apart from possible 
planets, a debris disk system contains remnant planetesimals 
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and dtist int o which they are ground through collisions (see, e.g. 
IWvatll2008l and references therein). Debris disks are expected to 
be nearly gas-free, at least extremely gas-poor compared to pro- 
toplanetary disks. Even so, i n the case of P Pic, gas was detecte d 
very early on in absorption (Slettebak 19751; iHobbs et al.lll985h . 
and later in emission (Olofsson et al. 2001), due to the favorable 
edge-on orientation of the disk. The observed gas around /3 Pic 
is most likely replenished, i.e. secondary, as opposed to a rem- 
nant from the initial star-fo rming cloud. Evidence for this comes 
from the presence of CO dVidal-Madiar et aljri994l; iJoUv et all 
ii998; Rober ge et al . ,20()0|). which would be dissociated on time- 
scales of -200 vr dvan Dishoeck & Blacklll988t iRoberge et all 
120001) and from the presence of ne utral gas elem ents in the 
disk dOlofsson et alJ 12001; Brandek eretal| |2004 . subject to 
short removal times (Ferna ndez et al. 2006). Possible mecha- 
nisms for producing secondary gas in clude photon-induced des- 
orption from solids ( Chen et ajl 20071) and grain-grain collisions 
dCzechowski & Mannll2007l) . Part of the observed gas may also 
stem from comet evaporatio n, as inferred from observed time- 
variable absorption lines (e.g. lFerlet et al.ll987l ; lBeust & ValironI 
,2007,) . 
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However, in general observations of gas are much more dif- 
ficuh than that of dust. Standard detection techniques either use 
CO as a tracer of hyd rogen which c an be observed at radio fre- 
quencies, as done bv iHughes et alj (|2008) for 49 Cet, or mea- 
sure H2 emission hnes which are pumped by stellar emission 
lines originating from the chromos pheric and coronal regions, 
which was done for AU Mic bv iFrance et alJ (12007 ). A poten- 
tially more sensitive way of finding gas is to look for it in ab- 
sorption, as was done for p Pic. The downside is that this re- 
quires the special edge-on geometry of the d isk, but th i s meth od 
has nevertheless been successfully used by iRedfieldl (|2007|) to 
detect circumstellar Nal absorption towards HD 32297, a star 
with a known disk. Conversely, stars which are known to ex- 
hibit circumstellar absorption lines, so called shell stars, can 
be searched for evidence of c ircumstellar material, as done by 
iRoberge & Weinbergeil (l2008l) using Spitzer/MIPS data. Out of 
16 surveyed shell stars they found infrared excess, and thus 
evidence for circumstellar dust, around four stars: HD 21620, 
HD 1 18232, HD 142926 and HD 158352. 

Despite substantial efforts, the gas component of the debris 
disks remains much less constrained observationally than the 
dusty one. It is quite possible that primordial gas survives longer 
than usually assumed, at least in the outer parts of the disks, or 
is present in larger amounts than expected, without showing up 
in observations. In fact, about ten Earth masses of gas, if not 
more, could still remain in many young debris disks where gas 
was searched for a nd not found, without violating observations 
( lHillenbrandll2008 b . If present, this hardly detectable gas would 
heavily affect the disk's physics and evolution and could neces- 
sitate revisions to standard theories of disk evolution and planet 
formation. 

The goal of our work is to analyze the effects of gas on the 
dynamical evolution of solids. We would like to find out whether 
gas, if present in larger amounts or with a different radial distri- 
bution than usually assumed, would alter the dust distribution 
and thus the brightness profile of a debris disk in such a way 
as to show u p in the observations. We foll ow the approach first 
suggested bv lThebault & Augereaul (l2005h who applied it to the 
the p Pic system: we first postulate a certain amount and spa- 
tial distribution of gas in one or another debris disk system, then 
compute a steady-state distribution of dust in it, calculate the 
observables such as brightness profile, and compare them with 
available observations. 

In Sect. 2 we select and analyze three young debris disk sys- 
tems relevant for this study. Sect. 3 lays down basic theory of the 
dust production and dynamical evolution in a debris disk with 
a gas component. Sect. 4 describes numerical simulations and 
Sect. 5 their results. In Sect. 6 we devise an analytical model 
and use it to interpret the numerical results. Sect. 7 contains our 
conclusions. 



2. Systems 

2.1. Selection criteria 

We wish to choose several young debris disks in which the pres- 
ence of gas in little to moderate amounts has been reported. 
Ideally, these should be edge-on systems, so that better con- 
straints on gas are available from the presence or absence of 
absorption lines, not just CO mm emission. We need dust disks 
that are spatially resolved, preferably in scattered light, so that 
the radial profile of brightness is known. The age of the disks 
should not be very far from the boundary that separates gas- 
and dust-rich, optically thick protoplanetary disks from nearly 



gasless, optically thin debris disks, which is believed to lie at 
a; 10 Myr The best ages would thus be 10-30 Myr 

We find three systems to satisfy these criteria the best: p Pic, 
HD 32297, and AU Mic. Known facts and key parameters of 
these systems relevant to our study are presented in the subse- 
quent sections. We stress, however, that all three systems should 
be regarded as "typical" examples of their classes and might be 
used as a proxy for other systems. Thus our results, being of in- 
terest in the context of particular objects, could at the same time 
be considered as generic. 

2.2. pPic 

Star. A 12 Myr-old (IZuckerman et al.ll200ll) A5V star at = 
19.44 + 0.05 pc. 

Dust and parent bodie s. Th e debris disk was first re- 
solved by Smith & Terrile (119841) and later at various wave- 
le ngths (A rtvmowicz 2000, and re ferences therein). According 
to I Artvm owicz & Clampin (Il997l) . the vertical optical depth 
of the dust disk has a maximum of 1.53 x 1 0'^ at 60 AU, 
and h as a slope of -1.7. ..-2 in the outer part. iMouillet et al.l 
d 19971) give 5 x 10 at lOOAU with an outer slope of 
-1.7. Dust mass is roughly 0.05...0.5Me dArtvmowiczl (20001: 
i Thebault & Augereaul l2005|). wit h O.l Mm bei n g probably the 
best estimate (Zuckerman & Becklin 1993; La grange et alJ 
200Q). Analysis by Augereau et al. (2001) (their Fig. 1) show 
that an extended dust disk produced by the planetesimal belt 
("birth ring") between 80-120 AU would closely match the re- 
solved scattered light images together with the long-wavelength 
photometric data. The dust distribution itself is given in Fig. 2 of 
lAugereau et al.l (l200lh . The radial SB profile of the mid-plane 
scattered -light images shows a slope of -3... -4 outside 130- 
260 AU (iGolimowski et al.l l2006'). The dust disk around p Pic is 
famous for its large scale asymmetries, which might be caused 
by a sub- stellar companion in the disk (see, e.g. [Mouillet et al.l 
ll997HAugereau et al.ll200ll) . 

Gas. Before ySPic was k nown to harbor a debris disk, it was 
classified as a shell star by 'SlettebalJ (1 19751) due to its promi- 
nent Ca II H & K absorpt i on. G as was then re-discovered in ab- 
sorption by iHobbs et al.l ( 1985|) an d in spatially resolved Na I 
emission by ' Olofsson et al.l (12001 1) . The atomi c hydrogen c on- 
tent of the disk was constrained by iFreudling et al.l (1 19951) to 
be < 2 Me, and the molecular column density to be < 3 x 
10'** cm"^ by iLecavelier des Etangs et al.l (i20()]J, which corre- 
spon ds to < 0.2M{B, assurn i ng the gas to be d istribu ted in the 
disk (iBrandeker et all l2004 . Brandeker et al.1 (|2004 observed 
spatially extended gas emission from a number of elements (in- 
cluding Na I, Fe I, and Ca II), and derived a spatial distribution 
for the gas 
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with «o - 2.25 x lO^* cm"^ for an assumed solar composition, and 
ro - 117 AU, leading to a total gas mass of 0.1 M^. They also 
investigated a metal-depleted case no = 10'' cm"-' (implying a 
total gas mass of 40 M®), which they found to be in contradiction 
with the observed upper limits on hydrogen. 



2.3. HD 32297 

Star. A 30 Myr-old A5V star (iManess et al.1 l2008l) at d = 
113 + 12 pc. 
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Dust and parent bodies. The dust disk was first resolved 
with HST/NICMOS in scattered light bv ' Schneider et all (l2005h 
up to 400 AU. The surface brightness (SB) of the SW wing 
is fitted by a power law with index -3.6, while the NE side 
shows a break at 200 AU: t he inne r part has a slope of -3.7, 
whereas the outer one -2.7. iKalasI (|2005) resolved the disk in 
the R-band between 560 and 1680 AU. The mid-plane slopes 
were found to be -2.7 and -3. 1 for NE and SW wing s, respec- 
tively, with strong asvmmetries. lMoerchen et an ( l2007l) resolved 
the disk with Gemini South/T-ReCS in therma l emission at 12 
and 18 j um up to 150 AU. Resolved images bv iFitzgerald et all 
(2007b) taken with Gemini North at 11 jim revealed a bilobed 
structure with peaks at ~ 65 AU from the star Maness et al.l 
(|2008|) marginally resolved the disk with CARMA at 1.3 mm. 

The spe ctral energy distribution (SED) fitting by 
IFitzgerald et al. (2007b) suggests a population of larger grains - 
and therefore a location of the birth ring - at w 70-80 AU. 
The vertical o ptical depth of the d ust disk at the same distance 
is 4 X 10"^ (iManess et al] l2008l) . 1.3 mm measurements by 
iManess et al.l (1200 8') point to the existence of a third population 
of even larger grains at a characteristic stellar distance of 50 AU, 
which probably comprises > 95 % of the total dust mass. Dust 
mass required to fit the SED up to far-infrared wavelengths is 
roughly 0.02 Mm, but 1 .3 m m flux may require as much as 1 
of dust (Man ess eta"ni2008l) . 

Gas. Redfieldl (l2007h found an intriguingly strong Na I ab- 
sorption. Assuming the morphology and abundances of the sta- 
ble gas component to be the same as for p Pic, and that the gas 
disk extends up to 1680 AU as debris disk does, he derived the 
total gas mass of ~ 0.3 M^. The absence of the observable CO 
/ = 2 - 1 emission with C ARMA places an upper limit on the 
gas mass of ~ 100 (Maness et al.ll2008l) . 



2.4. AU Mic 

Star. A 12 Myr-old dMle flare star, a member of the 
j3 Pictoris Moving Group ai d - 9.94 + 0.13 pc. It is the clos- 
est known debris disk resolved in scattered light. 

Dust a nd parent bodies. T he de bris disk was first resolved in 
R-band bv lKalas etall ( |2004 and lUu. (2004) . Later on, it was 
resolved with HS T/ACS by[ Krist et al,| j2005l) and in the H-band 
with Keck AO bv lMetc hev et al. (.20(351) . The dust fractional lu- 
minosity is 6 X 10""* |Liu 2004). Th e dust mass (up to 1 m m) 
is estimated to be ~ 2 x 10""*^/® ('Augereau & Beust 2006), 
but sub-mm fluxes require 1 x lO^^M® (Liu et al. 2004). The 
birth ring of planetesima ls is believed to be located at 35 AU 
dAugereau & Beustll2006l) . An R-band SB profile slope of -3.8 
between 3 5 -200 AU was fou nd by Kalas et al. (2004), whereas 
iLiul (120041) . Ikrist et al.1 (l2005h . and IFitzgerald et al.1 (I2007ah de- 
rived SB slopes in the range -3. 8. ..-4.7. Like jSPic and HR 
32297, the disk of AU Mic possesses asymmetries, which are 
probably formed by the dynamical influence of planets dLiul 
12004). 

Gas. Non-stringent upper limits on the gas mass wer e found 



from a search for optical absorption fines from Cal, Call and 
Fe I bv lBrandeker & Javawardhanal d2008l) . 



from non-detection of CO 3-2 emission by 'Liu et al.' ('2004') 
(< 1.3 Me) and H2 UV absorption by Roberge et al. (2005) 
(< 7 X 10'^ Mffi). [France et al.1 (12007 ) tentatively detected and 
analyzed fluorescent H2 emission. Within the observational un- 
certainties, the data are consistent with gas residing in the de- 
bris disk, although other possibilities such as a cloud that ex- 
tends beyond the disk cannot be completely ruled out. They 
found a very low total gas mass between ~ 4 x 10"'* M® and 
~ 6 X 10"^ Me, consistent with upper limits < 10""* Me obtained 



3. Basic theory 

3. 1 . General picture 

Throughout this paper, 
scenario of a debris 



we adopt the followi ng standard 



disk e volution (e.g. iKrivovetaLl 



2006t IStrubbe & Chiang! l 2006t iThebault & Augereaul 120071: 



Krivov et al 12008; Thebauh & Wull2008l) 



There is a relatively narrow belt of planetesimals ("birth 
ring") in orbits with moderate eccentricities and inclinations. 
We assume that this birth ring is located where the scattered 
image brightness peaks. Note that the systems resolved at 
(sub)-mm wavelengths usually exhibit a bright ring of ap- 
proximately the same radius. 

Orbiting planetesimals in the birth ring undergo collisional 
cascade that grinds the solids down to dust. We assume that 
the dust grains with radii [s, s + ds] are produced in the birth 
ring at a constant rate Nds, where 



(2) 



The parameter q is unknown. However, a usual assump- 
tion — which we will follow unless stated otherwise - is 
q = 3.5. 

At smallest dust sizes, stellar radiation pressure effectively 
reduces the mass of the central star and quickly (on the dy- 
namical timescale) sends the grains into more eccentric or- 
bits, with their pericenters still residing within the birth ring 
while the apocenters are located outside the ring. As a result, 
the dust disk spreads outward from the planetesimal belt. The 
smaller the grains, the more extended their "partial" disk. 
The dust grain orbits undergo slower modifications due to 
gas drag and experience gradual loss due to mutual colli- 
sions. 



3.2. Stellar gravity and radiation pressure 

We require that the disk is optically thin, so that each dust grain 
is fully exposed to stellar radiation at any location in the disk. 
Since the radiation pressure is proportional to r"^, as is the stellar 
gravity, a dust grain experiences "photogravity", i.e. gravity of a 
star with an "effective stellar mass" Mgff : 



Meff =M,(1-/?), 



(3) 



where /3 is the ratio of radiation pressure to gravity dBurns et al.l 
IT979I) : 



j3 = 0.5738 Qp, 



1 gem W l jLfm \ Li,/LQ 
Pbuik /\ s j M^,/Mq 



(4) 



Here, Qpi is the radiation pressure efficiency (henceforth set 
to unity), pbuik the material density of particles (we set it to 
3.3 gcm"^), s their radius, and and M* are stellar luminos- 
ity and mass, respectively. 

3.3. Gas drag 

We assume that the gas distribution remains unaffected by that of 
dust and that gas simply exerts a drag force on the dust particles. 
If the gas mass is larger than the dust mass, this assumption is 
natural. In the case where both are comparable, the validity of 
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this assumption will be checked later a posteriori. Indeed, we 
will choose the same initial distributions for dust and gas and 
will see that these will not diverge considerably in the course of 
the disk evolution. 

Gas orbits the star at a sub-keplerian speed 



(5) 



where v'k is the circular Keplerian velocity and 77 is the ratio of 
the force that supports gas against gravity to the gravity force. 

There exist two possible reasons for sub-keplerian rotation 
of a gas disk. One is the case of a thermally-supported disk, 
in which the sub-keplerian rotation stems from the gas pressure 
gradient. However, around /? Pic the gas at non-solar composi- 
tion has been observed, which is domin antly supported by ra- 
diation pressure rat her than gas pressure dFernandez et al.ll2006l: 
iRob erge et al. 2006). Thus we also consider another case, where 
gas is supported against stellar gravity by radiation pressure. 

For a thermally-supported gas disk, we follow a stan- 
dard description of the dust aerodynamics dWeidenschillingl 
Il977h. generalized to t he pre sence of radiation p r essure 
(Takeuchi & Artvmowicz 200?; iThebault & Augereaul 120051; 
[Herrmann & Krivov..2007.) . The factor 77 is 



Pg dr 



(6) 



Here, GM* is the gravitational parameter of the central star and 
P the gas pressure, proportional to gas density p„ and gas tem- 
perature Tg! 



■Pi 



jUgmn 



(7) 



with the Boltzmann constant k, the mean molecular weight fig 
and the mass of a hydrogen atom mn. 

The gas temperature and density are usually taken to be 
power laws 



and 

Pg r~^- 

Then, Eq. ^ takes the form 



(8) 



(9) 



(10) 



where is the gas temperature at a reference distance tq. The rj 
ratio only depends on the gas temperature slope p and its density 
slope ^, but not on the gas density at a given distance. (Note 
that sub-keplerianity requires p + 4 > 0.) As oc L^"*, the 

dependence on the stellar parameters is weak: rj oc L^J* /Mi,. 

We now consider the case where gas is supported against 
stellar gravity by radiation pressure. Denoting by jSgas the effec- 
tive radiation pressure coefficient acting on the gas, its speed is 



Vg = VK 



^1 -ySgas, 



yielding a simple relation 

=^gas- 



(11) 



(12) 



Regardless of the mechanism that supports the gas disk 
against gravity, the gas drag force on a dust grain is expressed 
by (Takeuchi & Artvmowicz 200 



Fd - -rrpgS^ (vj + Av^)''^ Av , 



(13) 



which combines the subsonic and supersonic regimes. Here, 
Av = Vd - Vg is the difference between the dust velocity Vd and 
the gas velocity Vg, and vj is the gas thermal velocity: 



4 / 8kTg ^'^2 



(14) 



3 yTTfiginu 

For later discussions of the timescales, we also define the stop 
ping time, 

4pd 5 1 



Stop 



(15) 



the time interval over which Av would be reduced by a factor of 
e if the drag force were constant. 

3.4. Gas temperature 

As we saw in the previous section, dust dynamics is expected 
to depend sensitively on the gas temperature, in particular on 
its radial gradient p. A commonly use d assumption is that the 
gas shares the dust temperature profile (iKamp & van Zadelhofll 
2001), which in the simple blackbody approximation gives p = 
1/2. This is a reasonable assumption if gas-dust interaction is 
strong and the photo-electric heating weak, but may not be 
valid in general. Indeed, in case of strong UV environments 
the photo-electric effect on dust can be the domin ant heating 
source of the gas and lea d to a dust drift instability ( Klahr & LinI 
l2005l:lBesla & Wull2007,) . For sufliciently high dust content, one 
may expect (at every distance) a power-law relation between 
the gas temperature and the number density of dust. Following 
Klahr & Lin ( 20051 Eq. (O generahzes to 



(16) 



with y > 0. A first-order heating-cooli ng analysis suggests y ~ I 
as a typical value (see the Appendix of lKlahr & Linl2005h . 

To evaluate how valid the p = 1/2 assumption is for a more 
detailed model of the thermal balance, we used the code Ontario 
(Zagorovsky, Brandeker, & Wu, in prep.). Ontario is tuned to 
model gas emission during the debris disk phase, given input pa- 
rameters related to the gas/dust disk structure, elemental abun- 
dances, and the stellar luminosity spectrum. Ontario computes 
the ionization and thermal balance self-consistently, with par- 
ticular care taken of heating/cooling mechanisms (the most im- 
portant being photo-electric and ionization heating, and cooling 
by C II 158 yum). The major simplifying assumptions are that the 
gas is considered to be in atomic/ion form (no molecules), and 
that the disks are optically thin (i.e. no chemistry and simplified 
radiative transfer), conditions that are expected to be closely met 
by debris disks around A and F stars, but not necessarily around 
stars of later spectral type. Using the same gas and dust pro- 
files as in the dynamical simulations (i.e., surface density slope 
of -1.5 corresponding to a mid-plane density slope of -2.5), we 
computed the mid-plane temperature for three different cases, as 
shown in Fig.[T](top): 

1 . an Mgas =0.1 Me disk at solar abundance with hydrogen 
entirely in atomic form (giving the mean molecular weight 
//g = 1.3) 
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2. a similar model but with Mgas - 10 

3. an Mgas =0.1 model with ySPic abundances (i.e. solar 
abundance except 20x carbon, no helium, and 10"^ x hydro- 
gen, giving /ig = 11.3, as motivate d by the inventory of gas 
observed around/? Pic, compiled bv lRoberge et al.ll2006l) . 

For comparison, the Tg oc r~'^^ is overplotted (with arbitrary 
normalization). The bottom panel of Fig. [1] shows the local p - 
-dlogT/dlogr exponent, which can be compared to the p - 
1/2 value. The numerical noise seen in the plot is due to the 
limited precision of the code {dT ~ 0.1 K). The temperature is 
not a power law, but varies smoothly with a local exponent that is 
constrained to the interval -1...0; a p = 1/2 power law therefore 
seems to be a reasonable approximation. 
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Fig. 1. Top: Gas temperature in the /3 Pic disk as a function of 
distance from the star, for three different assumptions about gas 
composition and density. (Dashed, dash-dotted, and solid lines 
correspond to tns, ths, and rns models, respectively, of Table [T] 
below.) Bottom: Radial slope p of the gas temperature. 



3.5. Radial dynamics of dust under piiotogravity and gas 
drag 

We start with a thermally-supported gas disk. Since gas adds 
or removes angular momentum to or from solid particles, it 
causes them to spiral outward or inward, until a certain size- 
dependent stability distance is reached, at which the gas pressure 
gradient and the stellar radiati on pressure balance each other 



graaient ana tne stellar raaiati or 
dTakeuchi & Artymowici200lh : 



;6(.9) = /7(i,rstab), (17) 

which can be solved for rstab for a given s or vice versa. 
Note that additional radial forces such as photophoresis 



dHerrmann & Krivovll2007l) can be included by simply adding 
them to the left part of Eq. ( [TT] ). 

We now turn to the case where the gas is radiation pressure- 
supported. Then, rj is independent of distance, so that there is 
only one value for rj over the entire disk and for all particle sizes. 
In this case, the particles with /3 < t] (J3 > rj) will spiral inward 
(outward) from the birth ring; those that just have P - rj will 
have no radial motion and will stay in the parent belt. Therefore, 
only grains with (3 > rj will make a contribution to the outer 
disk. They all will be drifting outward all the way through in a 
steady-state regime. 



3.6. Vertical dynamics of dust under piiotogravity and gas 
drag 

Apart from the radial dynamics described above, dust grain or- 
bits also have vertical evolution. The essential effect is dust sed- 
imentation (or settling) toward the mid-plane of the disk. It hap- 
pens on a timescale Tsett that depends on the particle size and the 
vertical profile of the gas density. Generally, Tsett also changes 
with time because of the grain's radial drift. If, for dust grains 
that contribute to visible brightness, the settling timescale is 
comparable to the stopping time (Tsett ~ 7"stop), a combination of 
settling and radial drift would cause the aspect ratio of the outer 
disk to decrease with distance ("anti-flaring"). Thus in principle, 
a comparison of the observed vertical distribution of brightness 
in an edge-on debris disk with a modeled one may offer another 
method of constraining the gas densities. In the present paper, 
however, we confine our analysis to the radial distribution and 
radial brightness profiles. 



3.7. Coiiisions 

Collisions are the main mechanism that limits the lifetime of dust 
grains in the outer disk, if they are large enough not to be blown 
away by radiation pressure or rapidly dragged away by gas. 

Collisional outcome is known to depend sensitively upon the 
relative velocities. Both proje ctiles are disrupted if their relative 
velocity Vrei exceeds (see, e.g. lKrivov et al.ll2005L their Eq. 5.2) 



2(mt + mp)2 

Vci- = ^ Q*D, 

V niinip 



(18) 



where nit and nip are masses of the two colliders and Q'^ is 
the critical energy for fragm entation and dispersal, which is 
~ 10^ erg g ' at dust sizes (e.g. lBenz & Asphaugll999 l). If Vrei < 
v'ci, the collision may result in partial fragmentation (cratering), 
restitution or, at low velocities, merging of the two colliders. 

The relative velocity Vrei mainly stems from the difference in 
the radial velocities of different-sized grains \vr\, as set by stellar 
photogravity and gas drag. In Sect. 5, we will see that Vrei is 
typically high enough for catastrophic collisions to occur even if 
the gas density is high. 



4. Numerical simulations 

4.1. Setups for simulations 

For each of the three systems, we adopt a fixed set of parameters 
for the central star and solids (planetesimals and dust) and test 
various gas density models (Table[T]i. One model may differ from 
another in four respects: 
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Table 1. Models for p Pic, HD 32297, and AU Mic. 
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* The arbitrary truncation distance (we set it to five times the central distance of the respective birth ring). It only affects the conversion of total gas 
mass to gas density. In that conversion, we also use a constant semi-opening angle of 15° for all disks, assuming the density to be constant with 
height. 



1. Gas support mechanism. We consider two possible rea- 
sons for sub-keplerian rotation of a gas disk: a thermally- 
supported disk, which seems reasonable in the case of pri- 
mordial gas at solar composition, and a radiation pressure- 
supported disk, which could be more app ropriate for sec- 
ondary gas with a non-solar composition. [Fernandez et al.l 
(I2OO6) computed the effective radiation pressure on the ion- 
ized gas to be y6 ~ 4 in the yS Pic disk (see their Fig. 4), 
assuming solar abundances. Since the mass is dominated by 
(inert) car bon, which was later found to be overabundant by a 
factor 20 (Roberge et al. 2006), the effective radiation pres- 
sure coefficient acting on the gas including 20x carbon is 
estimated to be ySgas ~ 4/20 - 0.2. Thus for the radiation 
pressure-supported disks we set 77 - fia-if, - 0.2. 

2. Radial slope of the gas temperature. In most of the models, 
we assume p - 1/2, but also test p - Q and p - \ Xo bracket 
the behavior of expected temperature profiles (Sect. [j4] i. 

3. Total amount of gas. In the nominal gas models, the total gas 
mass is taken as retrieved from the gas observations: O.IM® 
for p Pic, 0.3/V/e for HD 32297, 4 x 10 '♦/We for AU Mic. As 
we consider here only the outer part of the disks, outside the 
birth ring, we simply halve these masses. Indeed, in the yS Pic 
gas disk, the masses of the inner and outer disks are nearly 
equal (two terms in Eqs. (4) or (5) of Brand eker et al.l (|2004|) 
make comparable contributions). Interestingly, the above gas 
masses are roughly comparable with the dust mass (gas:dust 
ratio » 1:1). However, as explained in the introduction, the 
systems may contain much more gas than is evident in the 
observations. For this reason, we consider high gas mass 
models, in which the total gas mass is 100 times the nomi- 
nal mass. With this choice, the gas-to-dust ratio is a standard 
100:1. Finally, we try a very high gas case, in which the gas 
mass is ten times higher than in the high gas models. 

4. Radial slope of the gas surface density. If the gas is sec- 
ondary, one could expect t he gas profile to approxim ately 
follow the dust profile (e.g.. lCzechowski & Mannl2007l) . For 
our three systems, the latter falls off as oc r"'...oc slopes 
in this range are also expected on theoretical grounds in a 
standard "birth ring - collisionally evolving gas-free disk" 
model explained in Sect. 3.1. Thus our standard choice is to 



set the surface density of both gas and dust to be oc r"' ^ ini- 
tially. However, if the gas is primordial, i.e. a remnant of an 
accretion disk, the profile is more uncertain. One could still 
expect oc r"'-^ (consistent with an isothermal steady-state so- 
lution for a viscous accretion disk), but, for instance, a flat 
density profile oc r"' (another known steady-state solution) 
could also be possible. 

Accordingly, each model has an identifier XpYZ, where 

- X indicates the type of the gas disk 

- (no gas) 

- t (thermally-supported) 

- r (radiation pressure-supported) 

- p indicates the slope of gas temperature: 

-Oip^Q) 

- Hp^l) 

- nothing {p - 1/2) 

- Y indicates the total amount of gas: 

- (no gas) 

- n ("nominal" gas mass) 

- h ("high" gas mass) 

- V ("very high" gas mass) 

- Z indicates the slope of the gas density profile: 

- (no gas) 

- s (standard, surface density falls off as 1.5) 

- f (flat, surface density falls off as 1 .0) 

Our "tpns" and "rpns" models (a thermally- or radiation 
pressure-supported gas disk with various gas temperature pro- 
files, a nominal gas content, and a standard surface density slope) 
are expected to simulate secondary gas whose production is 
somehow related to dust. In contrast, the "ths", "tvs", and "thf", 
models would emulate a possible remnant of primordial gas in 
the outer part of the system. 

As reference models, we treat "tns" and "ths" cases (a 
thermally-supported gas disk with a nominal or high gas con- 
tent and a standard surface density slope). We ran these models 
for all three systems. In the case of /3 Pic, we additionally tested 
other disks, as listed in Tab. [1] 
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4.2. Numerical integrations 

We now describe the procedure to compute individual trajecto- 
ries of dust grains as well as the overall brightness profiles of the 
disks. 

First, we assumed that dust parent bodies are orbiting in the 
"birth ring". Their orbital semimajor axes were uniformly dis- 
tributed within the birth ring as specified in Table [T] and their 
eccentricities between [0.0,0.1]. 

In each run with a thermally-supported gas disk, we launched 
500 particles with radii distributed uniformly in a log scale be- 
tween s,„i„ and s,„a_^. The size ranges [s„,i„, s„,ax] used correspond 
to a j6-ratio interval from 0.9 down to 0.01 . The minimum /3 value 
of 0.01 is chosen to cover all grains which can have their stability 
distance in or outside the center of the birth ring (see Fig.|2]be- 
low). Upon release, the particles instantaneously acquire orbits 
with semimajor axes and eccentricities different from those of 
t heir parent bodie s, which is a standard radiation pressure effect 
dBurns et alJIT979l) . 

The forces included stellar gravity, direct radiation pressure, 
the Poynting-Robertson force, and gas drag. The drag force due 
to the stellar w ind, which is known to be important for AU Mic 
dStrubbe & Chi ang 2006; Augereau & Beust 2006), is not in- 
cluded. The reason is that a nearly radial stellar wind could no 
longer exist in the outer disk considered here, between roughly 
35 and 300 AU, because of the interactions with the presumed 
rotating gas. Indeed, a simple estimate, assuming a mass loss rate 
of ~ IO^'^Mq yr"' and the stellar wind velocity of ~ 400 km s"', 
yields the total mass of the stellar wind particles in the outer disk 
of $ lO^^Me, about three orders of magnitude less than the mass 

of rotating gas in the nominal model. 

The particle orbits were followed with the lEverhartl (11974 

Il985h integrator of 15th order with an adaptive step size. The 
integrations ended upon one of the following, whichever was the 
earliest: (i) a grain came as close as 10 AU to the star; (ii) a grain 
reached the distance of 1000 AU; (iii) after 10^ yr of integration. 
Instantaneous positions of particles were stored each 500 years 
for bound grains and each 5 years for unbound ones. A typical 
number of instantaneous positions per system and gas model was 
~ 10^ 

The setup for the runs with radiation pressure-supported gas 
disks was different from what is described above, because us- 
ing the same setup would lead to the following problem. As ex- 
plained in Sect. 4.3, we normalize the calculated dust density 
in such as way as to arrive at the correct maximum geometri- 
cal optical depth tq (Table [T]). In the usual runs for thermally- 
supported disks. To in the birth ring is dominated by the particles 
with j6 < 0.2. But in the "rns" run these drift inward, so that r 
comes from grains with /3 > 0.2. These have smaller cross sec- 
tion, and thus their number density turns out to be two orders 
of magnitude higher than in standard runs. Thus the collisional 
lifetime becomes quite short, ~ 1...10 years. For these reasons, 
the "rns", "rOns", and "rlns" models were run with just 15 grains 
(instead of 500) over 10^ years (instead of 10^), and the record- 
ing time step was as small as 0.5 years for bound grains and 
0.005 years for unbound ones. The minimum /3 value was set 
to 7/ = 0.2 (instead of 0.01), because in the radiation pressure- 
supported disks only grains with fi > rj drift outward from the 
birth ring. 

4.3. Collisional post-processing 

The collisions were applied to the numerical integration results 
through the following post-processing algorithm: 



7 

1. The instantaneous positions of grains stored during the nu- 
merical integrations are distributed into two-dimensional 
size-distance bins, (i,, rj) or (yS,, rj). The number of occur- 
rences in each bin Nunifii, rj) is converted into the absolute 
number density «(j6,-, rj) from the known normal geometrical 
optical depth at the birth ring (Table [T]l. Besides, for each 
bin we calculate the average radial velocity of its grains, 
Vr{l3i,rj). 

2. For each (/?,-, rj)-bm, the collisional lifetime of its particles is 
calculated as follows. We consider all bins with various ySt 
at the same distance rj and check if Eq.dTSll is fulfilled. The 
reciprocal of the collisional lifetime is 

Tcoiiifii,rj)-' = Y^n(J3k,rj) 

k 

X max[vv(j6/,rj),vv06t,r;)] 

X (T(J3i,l3k), (19) 

where cr is the collisional cross section and summation is 
performed over those projectiles that satisfy Eq.(fT8b. 

3. We then go back to the stored numerical integration results 
and go along each trajectory again. For a particle with y6, at 
each time step f/, we calculate a probability that the parti- 
cle survives destructive collisions over the current time step, 
p(J3i,ti) = e\p[-At/Tco\i(J3i,rj)], where rj is the distance at 
the moment f/ and At - tj- ti-\ is time step between two suc- 
cessive stored positions of the particle. Having p for all pre- 
vious time steps for the same grain, we determine the prob- 
ability for that grain to survive collisions up to the current 
moment of time: 

I 

P(J3i,ti)^Y[p^Pi^f'n^- (20) 

m=0 

4. Steps 1-3 are repeated iteratively. From now on, when 
determining the number densities in step 1, we use 
A^bin(y8;, rjfPiPi, t,) instead of just NbiM, rj). 

Note that we do not compute the exact relative velocities be- 
tween impacting grains, since this procedure would be too time 
consuming, but assume that this velocity is of the order of the 
vv of the grain having the highest radial velocity. In doing so, 
we implicitly neglect the azimuthal component of the relative 
velocity, and our procedure should thus be considered as giv- 
ing a lower estimate for impact speeds and thus shattering effi- 
ciency of collisions. The procedure is computationally fast and 
converges rapidly to a steady state. In practice, we perform five 
iterations. 

4.4. Calculation of surface brightness profiles 

From NbiniJSi, rj) recalculated after several iterations, the surface 
brightness of the disk is computed as 

SB{rj) = const x ^ A?bm(yS;, rj)s]-"rj' (21) 

(we set q - 3.5). Here, the factor s^^^ provides conversion from 
particle numbers to their cross section area (i^), accounts for a 
logarithmic binning of sizes {dN/dln s - sdN/ds), and includes 
the size distribution at production (i Eq.|2]i. In turn, the factor 
r"' is needed because the surface area of a radial annulus is pro- 
portional to r, whereas r"^ is a conversion from optical depth to 
\heSB. 
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The latter conversion needs to be explained in more detail. 
For the sake of simplicity, we assume here gray, isotropic scat- 
tering. In that case, the brightness of the edge-on disk is known 
to come mostly from the part of the line of sight c losest to the 
star (e.g. lNakanolll990HThebault & Augereaul2007l) . Further, as 
mentioned above, we assume a non-flared disk, with the dust 
scale height being linearly proportional to r. With these assump- 
tions, the standard brightness integral gives SB{r) oc T(r)r"^. 
Therefore, if 



T(r) oc r " 
then 

SB(r) oc r^' 
with 
- a + 2. 

5. Results 



(22) 



(23) 



(24) 



In this section, we start with an analysis of typical single-grain 
dynamics and then proceed with SB profiles. 

5. 1 . Radial drift of dust grains 

We start with the case of a thermally-supported gas disk. 
Numerical solution of Eq. (fTTT i with p - 1/2 is shown in Fig.|2] 
It shows that the stability distance is less than lOOOAU only 
for grains with /? ^ 0.15-0.20. Grains with higher/? ratios are 
sweeping outward all the way through the disk. 
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Fig. 2. Loci of stability in the /3 ratio - distance plane 
in thermally-supported gas disks (curves) and in a radiation 
pressure-supported disk (vertical straight line). Grains to the left 
of the stability lines drift inward, those to the right move out- 
ward. Thin horizontal lines mark the distances at which the birth 
rings around all three stars are centered. 

Fig. [3] illustrates how grains with different /3 ratios drift to- 
ward their stability distances in the nominal and high gas /3 Pic 
models. It shows in particular how sensitive the drift timescale 
2"drift (time needed for a grain to reach its stability distance) is to 
the the grain sizes and to the gas content. Further, it tells us that 
T^drift s> Tstop. One reason for this is that Tstop increases rapidly 
with increasing distance. For instance, in the nominal gas case, 
the y6 = 0.1 grains have the stability distance at ~ 500 AU, but 
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Fig. 3. Radial drift of different-size particles in the /3 Pic nominal 
gas (top) and high gas models (bottom). 



are still at ~ 150 AU after 1 Myr of evolution. Hence they have 
T'drift 3* 1 Myr, whereas their stopping time at the birth ring dis- 
tance is Tstop ~ 2 X 10^ yr. 

We now turn to the case where the gas is radiation pressure- 
supported. As explained in Sect. 13.51 only the particles with 
P > rj will spiral outward from the birth ring, and the stability 
distance does not exist. Thus the outer disk will only be com- 
posed of relatively small grains in unbound orbits. 



5.2. Radial velocities and collisional outcomes 

The relative velocity of impacting grains is a crucial parame- 
ter, leading to their destruction if it exceeds v^r (Eq. [TSl l. As 
discussed in Sect. 14.31 we assume that these impact speeds are 
of the order of particle radial velocities. Figure H] plots Vr for 
grains of various sizes, in the /3 Pic disk without gas and for two 
gas models, and compares them to the value of v„ ~ 300 ms"' 
for impacts between equal-sized particles. The top panel reveals 
Keplerian U-shape curves, where v,. > Vd for all grains smaller 
than s ~ 11 yum (J3 > 0.05), except near the apsides. Under the 
presence of gas (middle and bottom panels), radial velocities are 
damped, and only p ~ 0.1 {s ~ 10/im) grains have Vr exceeding 
the disruption threshold in the nominal gas case. For the high 
gas case, this damping is much more efficient so that after a few 
pseudo-orbits, even /3 ~ 0.2 (s ~ 4 yum) grains have < Vcr. The 
results are similar for the other two systems, i.e., low for all 
grains larger than p ~ 0.05 in the gas free case and p ~ 0.2 in 
the high gas case. 

However, collisional disruption plays a significant role in all 
the systems, for all gas models considered, and across the whole 
range of distances in the outer disk. This is because of the crucial 
role of collision between grains of unequal sizes. This is easy to 
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understand. Indeed, Fig. |4] shows that very small grains always 
maintain large radial velocities. Hence, for a grain of a given 
size, there is a significant amount of somewhat smaller grains 
that are fast enough and still massive enough to satisfy Eq. ( fTSl l. 
To illustrate this, we plot in Figure |5] the coUisional lifetimes of 
different-sized grains around all three stars, calculated with the 
algorithm described in section 5. We see that for /3 < 0.2, T^ow 
is more or less independent of particle sizes in all considered 
cases. Overall, Tcou is longer in high-gas models and shorter in 
nominal-gas ones. For bound grains (J3 < 0.5) in all models, they 
never exceed 4 x lO'* yr 



the one calculated with a full collisional simulation with ACE. 
A radial SB slope close to 3.5 is theoretically expected, too. The 
difference with the collisionless case is that the small high-/? 
grains can only be produced and destroyed in the birth ring but 
spend most of their time in collisionally inactive regions beyond 
it. As a result, their number density follows the ^ = 3.5 size dis- 
tribution within the birth ring (instead of being underabundant 
as in the collisionless case) so that the total integrated number 
of small grains (taking into account the large fraction outside the 
birth rin g) is much higher than the one derived from a q = 3.5 
law (see lStrubbe & ChiandllOOetlThebault & Wull2008h . 



5.3. Surface brightness profiles witlioutgas 

We now turn to SB profiles and start with a disk that does not 
contain gas ("000" model for Pic), both with and without col- 
lisions taken into a ccount. F or comparison, we ran our full col- 
lisional code ACE (iKrivov et al. 2008) that provides a detailed 
treatment of collisions with a multitude of collisional outcomes. 
The resulting profiles are shown in Fig.|6] 

Without collisions we get a steep -5) SB profile. This 
result is expected when assuming that grains down to the radia- 
tion pressure blow-out limit (J3 = 0.5) are produced in the birth 
ring with a q = 3.5 size distribution and that the smallest of 
these grains, which dominate the outer ring, are then simply di- 
luted along th eir eccentric orbits and th us get underabundant in 
the birth ring. lStrubbe & Chiang! (|2006|) showed that in this case 
the resulting surface density profile a pproximately decrease s as 
r -2 ^ yielding an SB slope of -4.5 (see lThebault & Wul2008l for 
a more detailed discussion). 

In contrast, the SB profile with collisions turned out to be 
close to -3. 5... -4.0, which corresponds to a surface density 
slope of X -1.5. ..-2.0. This slope agrees reasonably well with 
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Fig. 7. Same as Fig.|6] but for several gas models (thick curves). 
Collisions are included. 



5.4. Surface brightness profiles witli gas: p Pic 

Fig. |2]compares SB profiles for the yS Pic systems in several gas 
models. It shows that the nominal amount of gas leads to almost 
the same profile as the one without gas (close to -3. 5. ..-4.0). 
Furthermore, and surprisingly, all profiles with high and very 
high gas masses ("ths", "tvs", and "thf") are all close to each 
other and to those without or with little gas -3. 5. ..-4.0). 

Our "tvs" and "000" models are close, although not com- 
pletely identical, to the "hi gh gas" and "no gas" models in 
iThebault & Augereaul (l2005l) . In those two cases they obtained 
the slopes of ^ -2.5 and ^ -4.0, respectively (their Figs. 3 
and 4). Thus our result (w -3.0 and w -3.5. ..-4.0) slightly 
differs from theirs. Most of the difference comes from the fact 
that they used a more realistic, extended distribution of parent 
planetesimals, but a much simpler, monosized collisional model 
without a collisional disruption threshold. 

Next, we compared the "rns" profile with the "tns" one. Both 
turned out to be very similar (Fig. [8]l. This implies that both 
thermally-supported and radiation pressure-supported gas disks 
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Fig. 8. Same as Fig. |7] but for models assuming different gas 
temperature profiles. Since all the curves fall almost on top of 
each other, we artificially lowered each curve listed in the legend 
(starting from "tns") relative to the previous one by a factor of 
1.5. 



with the same amount of gas may yield similar radial distribu- 
tions of dust. 

Finally, we have also tested the influence of the radial pro- 
file of the gas temperature on the SB profiles. To this end, we 
compared the SB profiles of yS Pic obtained in the "tns" model 
{p - 111) with those in the "tOns" and "tins" models {p - 
and 1, respectively). Similarly, "rns" results were compared with 
"rOns" and "rlns". Again, the SB profiles turned out to be almost 
indistinguishable (Fig.|8). 

5.5. Surface briglitness profiles with gas: all systems 

We now proceed with the numerical runs for all three systems 
and two gas models ("tns" and "ths") for each system (TablelTJ. 
The resulting profiles are shown in Fig. |9] One lesson from the 
plots is about the role of collisions. In systems with a high dust 
density (or a large optical depth) and low gas density, collisions 
flatten the profile (J3 Pic, nominal gas; see also Fig |6] without 
gas). When the dust density is lower (AU Mic), collisions have 
little influence on the SB slope. The same is true in the case of 
a high gas density (high gas models for all three systems). This 
is mostly because the strong gas drag sustains sufficiently high 
radial velocities even far from the star, so that v,-(r) does not 
decrease with increasing r as abruptly as in the nominal gas cases 
(see Fig.m. 

In addition. Fig. |9] depicts partial contributions to SB pro- 
files by different-sized grains. The largest contribution typically 
comes from medium-sized grains with 0.05 < /3 < 0.5, most of 
which have stability distances outside the disk (Fig. |2| and can 
be treated as "effectively unbound"ones. The relative contribu- 
tion of the small /3 > 0.5 grains slightly rises with increasing 
distance, but never becomes comparable to that of the medium- 
sized particles. Large grains with /3 < 0.05 do not make any 
appreciable contribution to the SB profiles in any of the systems. 

However, the most important conclusion from Fig.|9]is that 
the slopes differ only moderately for the two extreme gas mod- 
els, an effect that we have already seen for /3 Pic and now see 
for the other two systems. In the nominal gas case, a distinc- 
tive feature of the profiles is their slight "curvature" — they do 
not follow a single power-law across the entire disk. This effect 
mostly comes from collisions rather than gas, and it can also be 
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Fig. 9. SB profiles for /? Pic (top), HD 32297 (middle), and AU Mic (bottom), for a nominal gas model (left) and a high gas one 
(right). Thick dashed and solid lines are without and with collisions taken into account, respectively. Their normalization is the same 
as in Fig. |6] Thin curves show partial contribution of three /^-intervals into the total profile with collisions. Note that shaded areas 
have another meaning than in Figs.|6]-[8] They now indicate an approximate range of profiles, as deduced from observations (see 
Sect. 2): -3.0...-4.0 foryS Pic, -2.7. ..-3.7 forHD 32297, and -3.8. ..-4.7 for AU Mic. 



seen in the gasless case (Fig.|6]l. The profiles are steeper close to 
the birth ring and are more gently sloping farther out. Between 2 
and 3rbirth (/"birth being the birth ring distance), the slopes are -4.6 
for /? Pic, -4.5 for HD 32297, and -4.7 for AU Mic. Outside 
3rbirth, the slopes flatten to -3.5, -3.2, and -3.3, respectively. In 
contrast, in the high gas case the curvature eff'ect is only present 
close to the birth ring. Outside 2rbirth, all three SB profiles have 
slopes in a -3.6. ..-4. Orange. 

Finally, the model slopes have to be compared with with the 
observed slopes: -3.0... -4.0 (/3 Pic), -2.7.. .-3.7 (HD 32297), 
and -3. 8. ..-4. 7 (AU Mic). From all these values and from Fig.|9] 
it is hardly possible to judge whether nominal-gas or high-gas 
models match observations better. In particular, this depends on 
the radial zone of the disk considered. Besides, one should keep 



in mind that our models rest on many simplifying assumptions 
(as, for instance, grey isotropic scattering) and have limited ac- 
curacy (e.g., contain some numerical noise). Equally, the slopes 
retrieved from observations inherit uncertainties from the data 
and are sensitive to the specific pro cedure of data reduction (see, 
e.g., a discussion in Sect. 3.2 of iFitzgerald et al] l2007al) and 
should be treated with caution, too. Thus the only conclusion 
we can make is that nominal-gas and high-gas models are both 
is reasonable agreement with observations. 



6. Analytic model 

To better understand the numerical simulation results, in this sec- 
tion we address the dust distributions analytically. 
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6.1. "Static" model in ttie case of a thermaiiy-supponed gas 
disk 

We start with a simple model that assumes dust to swiftly drift 
by gas drag: Tdiift Tufs, where Tdrift is the dust radial drift 
time and is the grain lifetime, e.g. collisional one. With this 
assumption, a dust grain of radius s is expected to "instanta- 
neously" arrive at an equilibrium distance ( fTTl l from the star r 
(for brevity the subscript "stab" will be omitted). This model 
ignores that fact that grains spend finite time on their way to sta- 
bility distances, and thus also contribute to the optical depth and 
brightness closer to the star than their parking distances. It also 
ignores that stability distances of smaller grains (if p < 1) are 
located outside the disk and they never arrive there. And, even if 
stability distances are well inside the disk, grains may not arrive 
there (or not all of them may), if they are collisionally eliminated 
on shorter timescales, meaning that the assumption Tdrift ^ T'life 
fails. 

Some of these assumptions fail in the systems considered in 
this paper. For instance, in Sect. 5.5 we showed that the largest 
contribution to the SB profiles comes from "effectively unbound" 
grains. Nevertheless, we deem this "static" model useful. First, 
the "static" case is fully tractable analytically. Second, it can give 
us a rough idea, to what extent the ambient gas can change the 
dust profiles and it can be considered as a limiting case for a 
more realistic "dynamic" model that will be worked out later. 

The grains with radii [s,s + ds] are produced in the birth ring 
at a constant rate Nds, see Eq. (|2]). These grains drift to a size- 
dependent stability distance, which is easy to find. In the limit of 
geometric optics, /3 oc Eq. Next, r] is given by Eq. ( fTOb . 
Equating /3 and rj we find a simple relation between the grain 
radius and stability distance (one-to-one if p 1): 



(25) 



Grains with radii [s, s + ds] are located in an annulus [r, r + dr]. 
Their steady-state number is 

dN = NTiif^ds, (26) 
and the normal geometrical optical depth of the annulus is 



T(r) 



dN 



Inrdr 



1 s(ry 



N{s{r))TM,(s(r)) 



(27) 



We now assume that the grain lifetime is independent of size (or 
distance). This is a reasonable approximation as can be seen, for 
instance, from Fig.|5] Substituting Eqs. ( |25] ) and (|2]) into Eq. (|27T i 
yields a radial dependence 



T(r) oc Tiife r " 
where 

a - 5 - 3p + (p - l)q. 
According to Eq. ( l24l i. the SB slope is 
- Ip + (p - \)q_ 



(28) 



(29) 



(30) 



In the above derivations we adopted the gas temperature 
which is independent of the dust distribution. However, as dis- 
cussed above (see Eq. [TSI l, in a high dust density limit the gas 
temperature may become proportional to a certain power of the 
dust density (or equivalently the optical depth t divided by dis- 
tance r). Using Eq. ( fTSI l instead of Eq. dS), Eq. (|29T l for the optical 
depth slope replaces by 



5-3p + {p-\)q-0-q)y 
\+y{3-q) 



(31) 



and Eq. ( [30] l for the SB slope by 

5-l>p + (p- \)q-0-q)y 



ju = 2 + 



\+yO-q) 



(32) 
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Fig. 10. Radial slope of surface brightness for different p, q, y, 
according to Eq. (|32] |. Thin and thick dotted lines (q = 3.0) co- 
incide with each other. Asterisk marks a standard model with a 
blackbody gas temperature {y = and p - 1/2) and a canonical 
size distribution of dust at production (q - 3.5). For compari- 
son, best-fit slopes inferred from the analysis of scattered light 
images of our three systems are shown on the right. 

Figure[TO]gives slopes p of surface brightness (without grain- 
grain collisions) for various p, q, and for y = and 1 . These re- 
sults can easily be understood. Consider, for instance, the stan- 
dard gas heating model (y - 0). We start with the dependence 
on q for a fixed p. If the temperature gradient is not too steep 
(p < 1), then s decreases with increasing r. In that case, assum- 
ing a steeper size distribution (larger q) makes the profile flatter, 
as it should. For a steep temperature drop-off (p > 1), it is the 
other way round. The steeper the size distribution, the steeper 
the profile. The dependence on p for a fixed q is also obvious. 
For all ^ > 3.0 (which is expected), a steeper temperature gra- 
dient steepens the SB profile. For all q < 3.0, the opposite is 
true. In the standard case that we took in most of the numerical 
simulations (p = 1/2, q - 3.5), we get p - 3.75 (asterisk in 

Fig. [13. 

In the case of a dust-induced gas heating (y - 1 ), Eq. ([32]iand 
Fig. [To] suggest that a much wider range of p is possible, and 
that the results are more sensitive to q. For example, the same 
standard case p = 1/2, q - 3.5 would result in an extremely 
steep p - 6.5. However, in view of the temperature calculation 
results presented in Fig.[T] we do not expect Eq. ( |32] | to be bet- 
ter approximation to reality than Eq. (i30] l. Rather, it is meant to 
show that the radial distribution of dust is quite sensitive to the 
assumed gas heating model. 



6.2. "Dynamic" model 

Instead of considering a system in which all grains reside at 
their stability distances, we now allow grains to drift through 
the disk towards their respective parking positions. We start with 
a thermally-supported gas disk and, for a time being, we ne- 



Krivov et al.: Gas in young debris disks 



13 



gleet eoUisions. Equation (|27J for the normal geometrieal optical 
depth of the disk replaces by 



T(r) 



- — ( ds — ( — 

2r Jo Vr(r, s) r Jo Vr 



-ds oc — I ds, 

r(r, s) r Jo vAr, s) 



and the surface brightness profile is now given by 



SB(r) 



- f 



Mr) 2-q 



Vr(r, s) 



-ds. 



(33) 



(34) 



Here, so(f) is the radius of those grains whose stability distance 
is r. For disks with p < I, sq decreases with r. The radius sq 
or equivalently, the radiation pressure-to-gravity ratio of grains 
with radius so, which we denote by /3o, can be read out from 
Fig.|2]or calculated analytically from Eq. ( fTOl l: 



/3o{r) = (p + 



(35) 



As a numerical example, for fi Pic at r = 500 AU we have Po ~ 
0.10. ..0.13 and sq ^ 6.6...8.6//m for all thermally-supported 
disks modeled. It is only grains with s S 6.6...8.6;um that are 
present at the distance of 500 AU. 

Unfortunately, y,.(r, s) is a complicated, non-power- 
law function of both argum ents (see, e.g., Eq. 23 in 
iTakeuchi & Artymowic3 l200lh . What is more, a distance- 
dependent integration limit so{r) in Eq. ( l34l l would result is a 
non-power-law brightness profile even if vv was a power law. 
The only straightforward particular case is blowout grains with 
yS t 0.5. Fig.|4]shows that their radial velocity is nearly constant, 
yielding SB oc r^^ . This is consistent with numerical simulations 
(see dotted lines in Fig.|9]). 

For the bound grains with moderate p that actually dominate 
outer disks, there are two competing effects. One is the distance- 
dependent integration limit so(r). Its physical meaning is as fol- 
lows. The larger the distance, the smaller the grains need to be 
to be able to reach it by gas drag. Thus larger grains are only 
present closer to the star, which steepens the SB slope. Another 
effect is that vv decreases with r (see Fig.|4|, making the SB slope 
flatter than -3. For instance, if we adopt v,- oc r"' as a very rough 
approximation for the outer disk, we will have SB oc r"^. 

Further complications are expected from collisions. The 
grains have then a limited lifetime, and the integrand in Eq. ( l34l l 
would have to be weighted with a fraction of particles that sur- 
vive collisions before they arrive at a distance r. This would gen- 
erally affect the SB slope. 

Because of the complexity of vv(r, s) and limited grain life- 
time due to collisions, it is difficult to extend our "dynamic" 
model further However, the model is still useful, as it uncovers 
the reason why SB slopes in systems with gas may be flatter than 
-3 (like in lThebault & Aug ereau 2005, see their Fig. 3): it is the 
slow-down of radial drift velocity with increasing distance. On 
the other hand, the model demonstrates that slopes steeper than 
-3 are also possible, because larger grains can only drift to lim- 
ited distances and only contribute to the parts of the disk close 
to the birth ring. 

Equations ([33]l-(l34li also hold for the radiation pressure- 
supported disks. However, in this case the upper integration limit 
So has another meaning: the radius of grains whose p ratio is 
equal to jSgas, Eq. ( fT2b . The above discussion applies in large 
part to this case, too. The main conclusion is that a slope around 
-3.0 is expected, possible deviations from which stem from a 
size-dependent radial drift velocity and collisions. 



7. Conclusions 

We considered young debris disks, in which there is observa- 
tional evidence for a rotating gas component (either primordial 
or secondary). We assumed that dust is replenished from parent 
bodies that are located in a "birth ring" (which usually shows 
up in the resolved images). We then modeled the dust distribu- 
tion and scattered-light brightness profile in the outer part of the 
disk, exterior to the birth ring, under different assumptions about 
possible amount and distribution of gas. 
Our main conclusions are as follows: 

1 . Our numerical simulations revealed that the radial profile 
of dust density, and thus the surface brightness profile of a dusty 
disk, are surprisingly insensitive to the parameters of a central 
star, location of the dust-producing planetesimal belt, dustiness 
of the disk and, most interestingly, the parameters of the ambi- 
ent gas. The radial brightness slopes in the outer disks are all 
typically in the range -3. ..-4. This result holds for gas densi- 
ties varying by three orders of magnitude and for different radial 
profiles of the gas temperature. Both the gas of solar composi- 
tion supported against gravity by gas pressure gradient and gas 
of strongly non-solar composition that must be supported by ra- 
diation pressure would lead to similar profiles, too. The slopes of 
-3... -4 are the same that were theoretically found for gas-free 
debris disks, and they are the same as actually retrieved from 
observations of many debris disks. 

2. Although the slopes roughly fall into the range -3. ..-4, 
the numerical simulations made it apparent that the exact slope 
depends on the total amount of gas in the disk and the gas den- 
sity distribution slope, as well as on the dust density (through 
collisional timescales). 

3. We developed a simple analytic description of the radial 
distribution of dust brightness in an optically thin disk. The an- 
alytic model explains the numerical results 1.-2. and provides 
guidelines to what can be expected in young debris disks, prior 
to any detailed numerical modeling. Assuming gray isotropic 
scattering by dust, the analytic model predicts a range of slopes 
around -3, due to the dominant contribution of small high-/? 
grains. It shows that deviations from this nominal value may 
come from the slow-down of radial drift of bigger grains at 
larger distances (flattening), from the fact that larger grains can- 
not reach larger distances (steepening), and from the collisional 
elimination of dust particles. In the limiting case of very high 
gas densities and low dust densities (where the "visible" dust 
drifts through the outer disk over timescales shorter than their 
collisional timescales), if the dust size distribution at production 
follows a power law with an index -3.5, and assuming a black- 
body gas temperature, the model predicts a slope of -3.75. 

4. Our results for three young (10-30 Myr old), spatially re- 
solved, edge-on debris disks {p Pic, HD 32297, and AU Mic) 
show that the observed radial profiles of the surface brightness 
do not pose any stringent constraints on the gas component of 
the disk. At least for gas densities falling within the observation- 
ally derived density limits, we do not see any significant effect 
of gas on dust distributions. Thus we cannot exclude that outer 
parts of the systems may have retained substantial amounts of 
primordial gas which is not evident in the gas observations (e.g. 
as much as 50 Earth masses for p Pic). However, the possibility 
that gas is only present in small to moderate amounts, as de- 
duced from gas detections (e.g. ~ 0.05 Earth masses in the p Pic 
disk or even less), remains open, too. In that case, gas would be 
secondary, stemming for instance from grain-grain collisions or 
photo-desorption of dust. 
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Appendix A: Nomenclature 

Symbol Description 



d 


Distance of the star from Earth 


G 


Universal gravitational constant 


k 


Boltzmann constant 


n 


Dust number density 


1 


Slope of dust size distribution in the birth ring (2) 


Ql 


Critical energy for fraginentation and dispersal )18t 


gp. 


Radiation pressure efficiency 


r 


Distance from the star 


t* 


Stellar luminosity 




Solar luininosity 


niH 


Mass of hydrogen atom 




Projectile mass jlS) 


nit 


Target inass )18t 




Stellar mass 


Mo 


Solar mass 


Me 


Earth mass 




Radius of a dust grain 


SB 


Surface brightness 


SED 


Spectral energy distribution 


P 


Radial slope of Tg 


P 


Gas pressure 




Gas temperature 




Collisional lifetime 119) 


Tdrift 


Time a grain needs to reach its rgtai, 


T'sett 


Dust settling time (towards the disk mid-plane) 


^stop 


Dust stopping time (15) 


Vcr 


Minimum relative velocity for collisional disruption ilS) 


Vg 


Orbital velocity of a gas parcel 


VK 


Circular Keplerian velocity 


Vr 


Radial velocity of a dust grain 


Vffil 


Relative velocity of two dust grains at collision 


Vt 


Thermal velocity of gas 


a 


Radial slope of t 


P 


Ratio of radiation pressure force to gravity for dust grains 




Ratio of radiation pressure force to gravity for gas 


y 


Efficiency of gas heating by dust (16) 




Ratio of gas-suppoi1ing force to gravity (6), ii2\ 


Pbulk 


Bulk density of dust grains 


Pg 


Gas density 


0" 


Collisional cross section of dust grains )19) 


T 


Normal geometrical optical depth 




Radial slope of SB 


Ms 


Molecular weight of gas 




Radial slope of gas density (9) 
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